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ABSTRACT:  Extensive  simulations  of  PAMAM  dendrimer  generation  2  were  performed  at  several  pH 
conditions  with  explicit  water  molecules,  to  obtain  proper  conditions  and  validity  for  additional  simulations 
without  explicit  water.  Within  the  range  of  validity,  simulation  without  water  greatly  extends  the  size 
and  duration  of  practical  simulations.  We  investigated  the  effects  of  long-range  interaction  parameters 
such  as  interaction  distance  and  dielectric  constant  for  molecular  dynamics  simulations  of  PAMAM 
dendrimer  without  water,  concluding  that  charged  dendrimer  simulation  with  distance-dependent 
dielectric  constant  but  without  cutoff  distance  best  mimics  explicit  water  results.  Structural  variations 
of  PAMAM  dendrimers  were  analyzed  as  a  function  of  pH  and  dendrimer  generation  using  MD  simulations 
with  these  long-range  interaction  parameters.  Globular  and  loosely  compact  structures  at  high  pH  (>  10) 
show  conservation  of  atom  density  distribution  across  dendrimer  generations.  Highly  ordered  extended 
structures  at  low'  pH  (<4)  present  an  increasingly  hollow  interior  as  dendrimer  generation  grow's,  resulting 
in  more  open  structure  which  provides  easier  access  by  chemical  agents.  By  contrast,  significant  branch 
back-folding  occurred  at  neutral  pH  in  addition  to  major  peripheral  distribution  of  the  terminal  groups, 
yielding  higher  interior  density  in  the  intermediate  radial  region  between  the  center  and  the  maximum 
radius  as  the  generation  grows.  Higher  generation  dendrimers  provide  a  cavity  surrounded  by  dense 
atom  populations,  producing  a  more  stable  agenL  carrier,  Transition  to  high-density  packing  occurs  between 
generations  4  and  5.  Volume  differences  between  neutral  and  low  pH  calculated  from  Rc  show  a  dramatic 
increase  beginning  at  generation  5. 


Introduction 

Synthetic  polymers  have  many  uses  in  biomedical 
applications,  based  on  their  ease  of  manufacture,  du¬ 
rability,  and  low  immunogenicity.1  Polyamidoamine 
(PAMAM)  dendrimers  are  the  only  polymers  with  uni- 
dispersed  and  well-defined  molecular  structures.  These 
molecules  can  be  synthesized  in  large  quantities  and 
have  a  large  number  of  potential  biomedical  applica¬ 
tions.  Several  groups,  including  our  own,  have  success¬ 
fully  used  PAMAM  dendrimers  as  a  gene  delivery 
system,2"4  and  we  are  now  developing  them  as  delivery 
systems  for  other  therapeutic  agents.  Other  groups  have 
found  PAMAM  dendrimers  suited  to  drug  delivery 
systems.50 

Desired  features  of  effective  drug  delivery  include  (1) 
protection  of  agents,  (2)  targeting  to  specific  tissues, 
cells,  and  cellular  compartments,  and  (3)  release  of 
agents  in  the  target  area  either  continuously  at  a  proper 
rate  or  in  response  to  external  signals.  There  are  several 
paths  of  uptake  of  exogeneous  materials  into  a  cell  with 
endocytosis  being  the  most  efficient  pathway  for  the 
uptake  of  macromolecules.7  Consequently,  escape  from 
the  endosome  is  critical  for  the  cytosolic  delivery  of  any 
macromolecule  before  it  is  degraded  in  the  lysosome. 
One  major  feature  of  the  endosome  is  its  local  pH  of  ~5, 
while  the  cytosolic  pH  is  ~7.4.  At  physiological  pH  7.4. 
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most  primary  amines  of  the  PAMAM  dendrimer  are 
protonated,  while  tertiary  amines  inside  the  dendrimer 
are  not  protonated  until  the  pH  drops  to  4.8  It  is 
speculated  that  changes  in  pH  structurally  alter 
PAMAM  dendrimers,  leading  to  endosomal  disruption. 

Computer  simulations  are  widely  used  in  polymer 
science.  Dendrimers  have  been  investigated  using  vari¬ 
ous  theoretical  models  and  computer  simulations  such 
as  self-consistent-field  calculations,9  Monte  Carlo  (MC) 
simulations,1011  and  molecular  dynamics  (MD)  simula¬ 
tions.1213  In  contrast  to  MC  simulation,  MD  simulations 
can  yield  dynamic  information.  There  are  two  general 
types  of  MD  simulations:  coarse-grain  and  atomistic. 
Coarse-grain  MD  simulations  can  give  information 
about  a  large  system  within  a  reasonable  computational 
time  frame,  while  atomistic  MD  simulation  offers  spe¬ 
cific  molecular  information  but  requires  considerably 
more  computational  time. 

While  coarse-grain  MD  studies  on  dendrimers13  were 
used  to  predict  general  properties,  atomistic  MD  simu¬ 
lations  on  dendrimers  were  performed  to  provide  3D 
structures  of  specific  dendrimers14"16  and  to  suggest 
possible  interactions  in  dendrimer  formation.17  Recently, 
more  exhaustive  atomistic  MD  simulations  on  dendrim¬ 
ers  were  carried  out  to  study  their  solution  behavior,18 
potential  use  as  catalyst  supports,19  the  effect  of  repeat 
unit  flexibility,20  and  several  nanoscale  applications.21 
Dendritic  polyelectrolytes  have  been  modeled  using  MC 
simulations,22  and  several  pH-dependent  dendrimer 
experiments  have  been  performed.23-24  However,  none 
of  these  simulations  have  dealt  with  pH-dependent 
protonated  dendrimer.  which  is  important  to  biological 
study  and  applications. 
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Figure  1.  Schematic  diagram  of  generation-2  dendrimer  at 
(a)  higli  pH,  (b)  neutral  pH,  and.  (c)  low  pH. 

Protonation  processes  at  different  pH,  a  part  of 
biological  events,  occur  in  aqueous  solution.  However, 
an  atomistic  MD  simulation  with  explicit  water  is  very 
expensive  in  both  computation  time  and  memory  and 
can  rapidly  exceed  the  capacity  of  current  software  and 
hardware.  Simulation  without  explicit  water,  on  the 
other  hand,  yields  uncertain  results.  There  are  methods 
to  mimic  a  water  environment  such  as  the  Poisson— 
Boltzmann  method,25  which  also  has  limitations  on 
model  size  and  takes  several  times  longer  than  simula¬ 
tion  in  vacuum  conditions. 

In  this  study,  we  used  atomistic  MD  simulations  to 
study  PAMAM  dendrimers  produced  from  an  ethylene- 
diamine  (EDA)  core  in  a  variety  of  simulated  physico¬ 
chemical  conditions  such  as  pH  and  dendrimer  genera¬ 
tion.  Since  simulations  with  water  molecules  are  very 
expensive,  we  investigated  parameters  for  simulations 
without  explicit  water  to  optimize  a  range  of  operation 
which  still  adequately  mimics  the  essential  character¬ 
istics  of  w'ater  effects.  Extensive  MD  simulations  of 
small  dendrimers  (generation  2)  were  performed  with 
explicit  water  molecules  at  different  pH  conditions. 
These  simulations  were  compared  to  simulations  with¬ 
out  water  while  changing  several  simulation  param¬ 
eters.  After  choosing  an  optimal  parameter  set  for 
simulations  without  water  molecules,  we  studied  pH 
effects  on  PAMAM  dendrimer  structures  as  a  function 
of  dendrimer  generation  under  these  optimized  condi¬ 
tions.  Structural  variance  corresponding  to  dendrimer 
generation  w'as  investigated  by  calculating  the  radial 
atom  density  of  the  dendrimer  from  the  center  of  mass 
as  well  as  the  terminal  nitrogen  atom  density  at  various 
pH  conditions.  In  addition,  the  structural  changes  of 
dendrimers  of  different  generations  in  low  (<4)  and 
neutral  pH  was  of  great  interest  because  of  pH  changes 
naturally  occurring  in  the  compartments  of  living  cells. 

Simulation  Conditions 


Table  1.  Physical  Properties  of  EDA  Core  PAMAM 
Dendrimer  Models 


generation 

pH* 

no.  of  atoms 

mol  wt 

total  charge 

2 

low 

546 

3286.5 

30 

neutral 

532 

3272.4 

16 

high 

516 

3256.3 

0 

3 

low 

1154 

6971.5 

62 

neutral 

1124 

6941.3 

32 

high 

1092 

6909.0 

0 

4 

low 

2370 

14341.6 

126 

neutral 

2308 

14279.1 

64 

high 

2244 

14214.6 

0 

5 

low 

4802 

29081.7 

254 

neutral 

4676 

28954.7 

128 

high 

4548 

28825.7 

0 

6 

low 

9666 

58561.9 

510 

neutral 

9412 

58305.8 

256 

high 

9156 

58047.8 

0 

a  Low  and  high  pH  refers  to  lower  than  pH  4  and  higher  than 
pH  10,  respectively.  Amine  protonation  depends  on  the  dendrimer 
generation  from  our  pH  titration  study  (data  not  shown). 


dendrimers  over  the  pH  range  to  ensure  the  hydrogen 
bond  interaction  range  was  covered.  After  the  steepest 
descent  minimization  process  for  5000  steps,  the  initial 
configuration  was  set  as  number  one.  MD  simulations 
were  performed  either  on  an  Onyx  using  a  consistent 
valence  force  field  (CVFF)2fi  of  Insight  II  or  on  a  Compaq 
AlphaServer  cluster  with  64  4-processor  nodes  using 
CHARMM27  software.  All  atoms  were  included  in  the 
potential  energy  calculations.  After  the  simulation 
continued  for  200  ps  at  295  K  (room  temperature)  with 
1  fs  intervals  from  the  initial  configuration  number  one 
(ICN1),  the  temperature  of  the  model  system  was 
increased  to  1000  K  and  equilibrated  for  1 1  ps  at  that 
temperature.  Among  the  configurations  of  the  1000  K 
simulation,  we  select  one  and  minimize  its  energy  for 
5000  steps  to  prepare  another  initial  configuration 
number  two  (ICN2).  Total  simulation  times  of  dendrim¬ 
ers  from  each  initial  configuration  were  400  ps  at 
neutral  and  low  pH  and  1.5  or  1.8  ns  at  high  pH.  Final 
configurations  from  ICNl  and  ICN2  were  compared  to 
confirm  equilibrium.  The  average  or  both  of  the  two 
configurations  were  used  for  comparison  to  results  from 
simulations  without  explicit  water  molecules. 

The  total  potential  energy  function  for  MD  calcula¬ 
tions  is  described  as 

^bonded  ^nonbonded 
^nonbonded  ^Lennard-Jones  ^Coulomb 


Based  on  pH  titration  data,  three  simulation  condi¬ 
tions  were  designed  to  represent  the  effect  of  pH:  no 
amines  are  assumed  to  be  protonated  above  pH  10,  all 
primary  amines  are  protonated  at  pH  7,  and  all  amines 
of  a  dendrimer  are  protonated  below  pH  4.  We  will  refer 
to  these  conditions  as  high  pH,  neutral  pH,  and  low  pH, 
respectively  (Figure  1).  Molecular  models  of  PAMAM 
dendrimers  (EDA  core)  from  generations  2  to  6  were 
built  for  these  pH  values  on  an  Onyx  workstation 
(Silicon  Graphics,  Inc.,  Mountain  View,  CA)  using  the 
Insight  II  software  package  (Accelrys  Inc.,  San  Diego. 
CA).  The  physical  properties  of  the  model  dendrimers 
are  shown  in  Table  1. 

Since  our  models  contain  strong  charges,  such  as 
those  found  in  aqueous  solutions,  it  is  necessary  to 
determine  the  validity  of  simulation  without  explicit 
water  molecules.  This  was  tested  using  dendrimers  of 
generation  2  with  a  10  A  water  layer  on  each  of  the 


(1) 


where  e  is  the  minimum  energy  of  Lennard-Jones 
potential,  u  is  the  distance  to  give  minimum  Lennard- 
Jones  potential,  q  is  the  partial  charge  on  the  atom,  D 
is  the  dielectric  constant  (1  for  vacuum),  ris  the  distance 
between  i  and  j,  and  1,  j  are  nonbonded  atom  pairs. 
Interactions  between  nonbonded  atoms  are  usually 
short  range  if  the  Coulomb  energy  term  is  not  dominant 
in  the  above  equation.  Therefore,  most  MD  simulations 
calculate  potential  energies  between  particles  only 
within  a  limited  nonbond  cutoff  distance  (r,).  However, 
since  our  models  at  low  and  neutral  pH  contain  many 
positive  charges,  and  the  long-range  Coulomb  energy 
is  dominant,  we  need  to  consider  the  effect  of  long-range 
interactions. 
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Figure  2.  Final  simulated  configurations  at  (a)  high,  (b) 
neutral,  and  (c)  low  pH,  starting  from  the  initial  configuration 
number  one  (left)  and  number  two  (right).  Water  molecules, 
explicitly  included  in  the  actual  simulation,  are  not  shown. 


Simulations  of  generation  2  (G2)  nonhydrated  PAM- 
AM  dendrimers  were  performed  under  the  two  com¬ 
bined  conditions  of  nonbonded  potential  energy  calcu¬ 
lation:  (1)  potential  energy  calculation  range  (nonbond 
cutoff  distance  rc)  and  (2)  dielectric  constant  D  in  the 
Coulomb  energy.  We  either  included  all  atoms  in  the 
energy  calculation  (rc  =  °°)  or  did  not  include  potential 
energies  between  atoms  farther  than  1 1  A  apart  (rc  = 
11  A).  We  also  used  the  dielectric  constant  D  —  1 
(vacuum)  or  distance-dependent  dielectric  constant  ( D 
=  r).  After  the  steepest  descent  minimization  for  1000 
steps,  each  model  was  annealed  for  5  ps  at  1000  K. 
Then,  the  temperature  was  cooled  to  295  K  and  simula¬ 
tion  performed  for  100  ps.  All  simulations  were  per¬ 
formed  on  an  Onyx  using  Insight  II  software  and  on  an 
AlphaServer  cluster  using  CHARMM  software  to  check 
software  and  hardware  dependency. 

On  the  basis  of  comparison  between  simulation 
results  with  and  without  water,  we  chose  the  long-range 
interaction  conditions  for  nonhydrated  PAMAM  den- 
drimer  simulations  such  as  (1)  rt  =  «  and  (2)  D  =  r. 
Under  these  conditions,  PAMAM  dendrimers  of  genera¬ 
tions  3-6  (EDA  core)  were  simulated  at  different  pH 
values  on  an  Onyx  using  the  CVFF  force  field  of  the 
Insight  II  software  package.  All  simulations  were 


distance  (A)  dlsUncs  (A) 

(<=)  m 


Figure  3.  Density  profiles  of  G2  dendrimer  with  water 
molecules  at  high  pH  (a,  b),  neutral  pH  (c,  d).  and  low  pH  (e, 
f).  Left  panels  (a,  c,  e)  show  the  atom  densities,  and  right 
panels  (b,  d,  f)  show'  the  densities  of  terminal  nitrogen  atoms. 
Blue  plots  are  results  from  ICN1  and  red  from  ICN2.  These 
graphs  correspond  to  Figure  2. 

performed  for  100  ps  (at  least  twice  the  equilibration 
time)  at  295  K  with  1  fs  intervals  after  energy  minimi¬ 
zation  and  system  annealing  processes  as  previously 
described.  To  ensure  the  equilibration,  the  larger  gen¬ 
eration  5  and  6  dendrimers  were  simulated  using 
another  initial  configuration  similar  to  the  G2  models 
with  water.  All  data  used  for  evaluation  were  collected 
from  the  simulations  after  equilibration. 

Results  and  Discussion 

The  computation  time  of  simulations  with  explicit 
water  molecules  is  extremely  expensive.  The  average 
service  unit  (SU:  1  SU  =  1  h  of  CPU  time)  for  100  ps 
simulation  of  G2  PAMAM  dendrimer  with  a  10  A 
hydration  layer  at  high  pH  and  1  fs  time  step  (ts)  was 
170.24  (5.32  h  with  32  processors)  on  a  Compaq  Al¬ 
phaServer  cluster  with  667  MHz  21264A  processors.  The 
SU  charge  for  the  100  ps  simulation  of  G2  dendrimer 
without  water  under  the  same  conditions  {D  —  1,  rc  = 
»,  ts  =  1  fs)  was  1.52  (0.38  h  with  four  processors), 
which  is  the  longest  among  nonhydrated  models.  There¬ 
fore,  simulations  with  water  are  more  than  100  times 
more  expensive  than  nonhydrated  simulations  for  den¬ 
drimer  generation  2.  This  difference  can  even  increase 
at  higher  generation  levels  limiting  the  current  practi¬ 
cally  of  explicit  water  simulations  of  dendrimers  for 
biological  applications.  If  we  can  predict  a  major  feature 
from  nonhydrated  simulations  with  a  proper  parameter 
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Figure  4.  The  G2  dendrimer  configuration  with  water  molecules  after  200  ps  simulation  from  ICN1  at  high  pH.  The  annealing 
process  at  1000  K  to  generate  ICN3  started  from  this  configuration.  Black  lines  represent  hydrogen  bonds  between  water  molecules 
and  the  G2  dendrimer  at  high  pH.  These  bonded  water  molecules  are  rendered  as  ball-and-stick,  while  all  other  water  molecules 
are  blue-lined. 


choice,  the  simulation  serves  as  a  useful  and  practical 
tool.  This  helps  to  guide  later  investigation  into  more 
specific  applications  that  may  require  explicit  water 
simulation  to  capture  specific  details. 

G2  PAM  AM  dendrimer  molecules  with  a  10  A  hydra¬ 
tion  layer  at  each  of  the  three  pH  conditions  were 
simulated  without  nonbond  cutoff  in  order  to  include 
all  long-range  charge  interactions  in  the  model.  ICNl’s 
at  each  of  the  pH  conditions  were  very  close  to  each 
other  (homogeneous  and  extended).  After  the  initial  200 
ps  simulation  and  consecutive  annealing  process,  we 
could  obtain  very  different  ICN2's  from  ICNl’s  at 
different  pH  conditions.  The  final  configurations  after 
MD  simulations  for  each  pH  are  shown  in  Figure  2. 
comparing  the  cases  of  ICN1  and  ICN2.  The  final 
configurations  of  neutral  and  low  pH  are  after  400  ps 
simulation  time.  Their  features  show  high  similarity 
between  two  different  initial  configuration  cases.  The 
final  configurations  of  high  pH  are  after  1500  ps 
simulations  yet  still  present  somewhat  different  features 
depending  on  the  initial  configuration.  Neutral-  and  low- 
pH  dendrimer  molecules  reached  equilibration  much 
faster  than  high-pH  molecules,  probably  due  to  exten¬ 
sive  charge  interactions.  In  the  high-pH  model  there  is 
little  driving  force  to  overcome  the  energy  barrier 
between  the  two  configurations  in  Figure  2a,  where  the 
net  charge  is  zero. 


We  calculated  the  atom  density  of  dendrimer  en¬ 
sembles  collected  from  snapshots,  after  system  equili¬ 
bration  was  assumed.  For  each  snapshot  configuration 
we  calculated  the  histogram  of  the  distances  of  all  atoms 
from  the  center  of  mass.  All  collected  histograms  of 
dendrimer  ensembles  were  normalized  by  the  total 
number  of  atoms  (unless  described  separately).  The 
density  of  terminal  nitrogen  atoms  in  amine  groups  was 
similarly  calculated  on  the  basis  of  radial  distance  of 
terminal  nitrogen  atoms  from  the  center  of  mass 
normalized  by  total  number  of  terminal  nitrogen  atoms 
(unless  described  separately). 

Figure  3  shows  the  atom  density  and  the  density  of 
terminal  nitrogen  atoms  corresponding  to  Figure  2. 
Each  graph  consists  of  500-1000  ensembles  close  to  the 
final  configuration.  The  final  configurations  at  low  pH 
give  a  very  narrow  atom  distribution  including  the 
distribution  of  terminal  nitrogen  atoms  localized  in 
finite  range  (Figure  3e,f).  There  is  a  slight  shift  of 
density  profiles  at  neutral  pH  (Figure  3c.d).  However, 
the  inner  radial  profiles  of  the  atom  density  are  very 
similar,  and  the  major  peak  of  terminal  nitrogen  density 
profile  of  each  model  is  found  at  the  same  place. 
Therefore,  we  can  assume  that  simulations  at  neutral- 
and  low-pH  conditions  are  equilibrated.  The  average 
values  from  the  two  initial  configurations  for  neutral 
and  low  pH  will  be  used  as  references. 
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(a) 


Figure  5.  Atom  densities  of  G2  dendrimer  with  water 
molecules  at  high  pH  as  time  progresses.  Each  graph  consists 
of  300-500  ensembles.  Simulation  results  at  0.3  ns  (blue),  0.6 
ns  (purple).  0.9  as  (green),  1.2  ns  (orange),  1.5  ns  (red),  and 
1.8  ns  (pink)  from  (a)  ICNl  and  (b)  ICN2. 

At  high  pH.  however,  density  profiles  from  ICN2  are 
more  compact  than  density  profiles  from  ICN1  (Figure 
3a, b).  Moreover,  the  general  shapes  of  the  profiles  are 
different  from  each  other.  Since  the  total  charge  of  high 
pH  models  is  zero,  the  hydrogen  bond  energy  (included 
in  Coulomb  energy  term  in  eq  1)  and  the  Lennard-  Jones 
energy  are  dominant  among  the  nonbonded  potential 
terms  in  eq  1.  Figure  4  shows  the  hydrogen  bond 
interaction  between  G2  dendrimer  and  water  molecules 
after  200  ps  simulation  from  ICN1  at  high  pH.  This  is 
just  before  the  simulation  at  1000  K  for  the  preparation 
of  ICN2.  During  the  simulation  at  1000  K,  hydrogen 
bonding  between  water  molecules  is  disrupted,  resulting 
in  the  water  being  dispersed  into  vacuum  space.  During 
this  brief  removal  of  water  the  major  energy  contribu¬ 
tion  comes  from  van  der  Waals  attraction  so  that  some 
portions  of  the  dendrimer  chains  are  able  to  rearrange. 
Therefore,  ICN2  is  compact,  while  ICN1  is  extended. 

We  tracked  the  simulation  progress  of  hydrated  G2 
dendrimer  at  high  pH.  Figure  5  represents  the  atom 
densities  of  G2  dendrimer  from  ICNl  and  ICN2  at 
progressing  time  intervals.  Each  plot  consists  of  300- 
500  ensembles.  The  ICN2  plots  stay  within  a  confined 
radial  distance  range,  indicating  a  local  minimum  from 
which  the  system  does  not  escape  within  the  simulation 
time.  The  ICNl  plots  move  to  the  left  on  the  distance 
axis  as  time  progresses  and  generally  settle  into  a 
different  local  minimum  except  in  the  outer  tail  portion. 
We  will  use  the  results  from  both  ICNl  and  ICN2  for 
the  high-pH  case,  understanding  their  differences  within 
our  simulation  limitations. 

Nonhydra  ted  C2  dendrimers  were  simulated  under 
several  long-range  interaction  conditions.  The  effect  of 
nonbond  cutoff  distance  r,.  was  of  major  interest.  Simu¬ 
lations  with  rc  =  11  A  and  rc  =  °°  were  compared.  We 
also  hypothesized  that  the  most  important  role  of  water 


Table  2.  Rc  of  PAMAM  Dendrimer  Generation  2 


water 

/; 

D 

PH 

Rci A) 

yes 

no 

1 

low 

18.87  ±0.11a 

neutral 

16.47  ±0.61'> 

high 

13.28  ±0.12,  9.47  ±  0.146 

no 

ll  A 

1 

low 

16.35  ±0.11 

neutral 

14.03  ±0.05 

high 

9.26  ±0.14 

no 

nO 

1 

low 

20.06  ±  0.09 

neutral 

17.44  ±0.11 

high 

9.16  ±0.18 

no 

ll  A 

V 

low 

15.30  ±0.18 

neutral 

13.05  ±0.21 

high 

8.28  ±  0.07 

no 

no 

V 

low 

16.59  ±0.13 

neutral 

14.47  ±0.34 

high 

8.43  ±0.10 

a  Average  value  between  results  from  ICNl  and  ICN2.  6 Each 
value  represents  results  from  ICNl  and  ICN2,  respectively. 


in  our  model  is  to  shield  charges  because  of  its  high 
dipole  moment.  We  changed  the  vacuum  dielectric 
constant  {D  =  1)  to  the  distant-dependent  dielectric 
constant  (D  =  r)  and  compared  the  results  with  the  data 
from  simulations  with  water. 

The  radius  of  gyration  Rc  of  all  G2  dendrimer  models 
was  calculated  as 


and  is  presented  in  Table  2.  All  models  show  the  same 
trends  in  pH  dependency;  with  lower  pH,  Rc  increases. 
The  difference  between  high  and  neutral  pH  is  larger 
than  between  neutral  and  low  pH.  At  neutral  and  low 
pH,  Rc  of  the  D=  l,  rc  =  1 1  A  model  is  similar  to  the 
D  =  r,  rc  =  °°  model,  the  latter  being  slightly  larger. 
Simulations  with  water  molecules  have  the  same  long- 
range  parameters  of  the  D—  1 ,  rc  =  °°  model.  Including 
explicit  water  reduced  the  size  of  Rq.  At  high  pH,  D  = 
1  models  show  values  closer  to  the  simulation  results 
of  hydrated  models  than  D  =  r  models. 

Figures  6,  7,  and  8  show  the  atom  densities  and  the 
densities  of  terminal  nitrogen  atoms  in  G2  dendrimers 
without  water  molecules,  compared  to  the  results  of 
simulation  with  water  at  high,  neutral,  and  low  pH.  For 
relative  comparisons  between  the  total  atom  and  the 
terminal  nitrogen  atom  densities,  the  densities  of 
terminal  nitrogen  atoms  were  normalized  by  the  total 
atom  numbers  of  their  models,  and  the  density  values 
of  terminal  nitrogen  atoms  were  increased  to  3  times 
their  real  values  for  illustrative  purposes.  Since  the 
terminal  amines  are  reaction  sites,  special  attention 
should  be  paid  to  the  distribution  of  terminal  nitrogen 
density  in  relation  to  the  total  atom  density  profile. 

Figure  6  shows  the  density  profile  from  simulations 
at  high  pH.  Simulations  from  ICNl  (Figure  6e)  have  a 
radius  more  similar  to  all  the  nonhydrated  models  than 
to  the  simulations  from  ICN2  (Figure  6f)-  The  atom 
density  profile  of  the  ICNl  model  parallels  the  increase 
of  the  D=  1 ,  rc  =  «  model  up  to  approximately  8  A  and 
then  from  12.5  A  decreases  closely  with  that  of  the  D  = 
1,  rc  =  1 1  A  model  (Figure  6g).  The  atom  density  profile 
from  the  ICN2  model  increases  closely  with  that  of  the 
D  —  r  models  up  to  5  A  but  then  diverges  from  the  other 
models.  None  of  models  without  water  show  a  peak 
shape  of  the  atom  density  plot  similar  to  the  ICN 1  or 
ICN2  models.  The  density  profiles  of  the  terminal 
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Figure  6.  Atom  densities  (blue)  and  the  densities  of  terminal 
nitrogen  atoms  (red)  in  G2  dendrimers  at  high  pH.  The 
nonbonded  interaction  conditions  of  models  without  explicit 
water  molecules  are  (a)  D  =  1  and  n  =  11  A.  (b)  D—  1  and  rc 
=  (c)  D  =  r  and  r,  =  11  A,  and  (d)  D  —  r  and  r,.  =  <*>.  They 

are  compared  with  the  results  of  simulation  with  water  from 
(e)  ICN1  and  (f)  ICN2.  To  improve  the  visual  comparison,  all 
density  profiles  of  terminal  nitrogen  atoms  are  increased  to  3 
times  their  real  value.  The  atom  densities  of  all  models  are 
compared  in  (g):  D—  1  and  r,:  —  11  A  (blue  line),  D—  1  and  rc 
=  oo  (pink  line),  D=  rand  n  —  11  A  (green  line).  D  —  rand  r,; 
=  «■  (red  line),  ICN1  (purple  circle),  and  ICN2  model  (brown 
circle). 

nitrogen  atoms  present  several  notable  features.  The 
ICN1  model  has  several  peaks.  The  D  —  1.  rc  =  11  A 
and  D  =  1.  rc  =  °°  models  (Figure  6a, d)  show  similar 
peak  distribution.  Based  on  these  density  profiles,  the 
distribution  range  of  the  ICN 1  model  is  between  the  D 
—  1 ,  rc  =  1 1  A  model  and  the  D  =  1 ,  rc  =  °°  model.  None 
of  the  nonhydrated  models  can  reproduce  the  ICN2 
model,  while  the  D  —  rc  =  1 1  A  model  is  the  closest 
mimic  of  the  ICN  2  model. 

The  density  profiles  of  simulations  at  neutral  pH  are 
presented  in  Figure  7.  From  the  D  =  1,  rc  —  °°  model 
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Figure  7.  Atom  densities  (blue)  and  the  densities  of  terminal 
nitrogen  atoms  (red)  in  G2  dendrimers  at  neutral  pH.  The 
nonbonded  interaction  conditions  of  models  without  explicit 
water  molecules  are  (a)  D  —  1  and  /;•  =  11  A,  (b)  D—  1  and  rc 
=  ■*>,  (c)  D=  r  and  rr  =  1 1  A,  and  (d)  D=  r  and  r,  —  They 
are  compared  with  (e)  the  results  of  simulation  with  water, 
which  are  the  average  values  of  results  from  ICNl  and  ICN2. 
To  improve  the  visual  comparison,  all  density  profiles  of 
terminal  nitrogen  atoms  are  increased  to  3  times  their  real 
value.  The  atom  densities  of  all  models  are  compared  in  (f) : 
D  —  1  and  rc  =  1 1  A  (blue  line) ,  D  =  1  and  r  —  °°  (pink  line) , 
77=  rand  /;:  =  11  A  (green  line),  D=  rand  rc  =  ~  model  (red 
line),  and  models  with  water  molecules  (purple  circle). 


(Figure  7b),  a  10  A  water  layer  changed  the  density 
distribution  dramatically  (Figure  7e).  The  atom  density 
distribution  (blue)  is  expanded  to  a  greater  radius,  and 
the  density  of  terminal  nitrogen  atoms  (red)  spread 
down  toward  the  center  of  mass.  An  additional  outside 
tail  bump  also  appeared,  which  is  not  present  in  any  of 
the  nonhydrated  models.  The  atom  density  profile  of  the 
with-water  model  increased  closely  with  the  D  =  r,  rc  = 
»  model  up  to  ~  12  A  and  then  spread  out  to  larger 
radius  (Figure  7f).  Nevertheless,  except  in  the  tail  region 
of  the  hydrated  model,  the  shapes  of  the  two  models  are 
similar.  Additionally,  some  terminal  nitrogen  atoms  of 
the  Z7  =  r.  rc  =  model  begin  close  to  the  center  of  mass 
and  a  main  peak  of  the  density  profile  of  terminal 
nitrogen  arises  close  to  the  outer  part  of  the  model 
(Figure  7d),  which  presents  a  very  similar  profile  to 
results  from  the  hydrated  model  (Figure  7e). 

Figure  8  shows  the  density  profile  of  simulations  at 
low  pH.  With  the  addition  of  alO  A  hydration  layer, 
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Figure  8.  Atom  densities  (blue)  and  the  densities  of  terminal 
nitrogen  atoms  (red)  in  G2  dendrimers  at  low  pH.  The 
nonbonded  interaction  conditions  of  models  without  explicit 
water  molecules  are  (a)  D  =  1  and  rc  =  1 1  A.  (b)  D  =  1  and  ry 
=  (c)  D  —  r  and  rc—  11  A.  and  (d )  D—  r  and  rc  =  «.  They 

are  compared  with  (e)  the  results  of  simulation  with  water, 
which  are  the  average  values  of  results  from  ICN1  and  ICN2. 
To  improve  the  visual  comparison,  all  density  profiles  of 
terminal  nitrogen  atoms  are  increased  to  3  times  their  real 
value.  The  atom  densities  of  all  models  are  compared  in  (f): 
0=1  and  rt  —  1 1  A  (blue  line),  D~  1  and  rt  =  ~  (pitik  line), 
D  =  rand  rc  —  11  A  (green  line),  D  =  rand  r,  =  °°  model  (red 
line),  and  models  with  water  molecules  (purple  circle). 

some  shrinking  of  the  radial  distance  from  the  D  =  1 , 
rc  =  °°  model  is  observed  (Figure  8f).  The  atom  density 
profiles  of  all  models  increase  similarly  up  to  12  A.  After 
12  A,  the  density  profile  of  the  hydrated  model  lies  in 
between  the  D  =  r,  rc  =  °°  model  and  the  D  =  1 ,  rc  =  °° 
model,  though  closer  to  the  latter.  The  largest  4  A  region 
of  atom  density  profiles  of  the  D=  r,  rc  —  °°  model  and 
hydrated  model  are,  however,  very  similar.  Interest¬ 
ingly,  the  two-peak  profile  of  the  terminal  nitrogen 
density  of  the  D=  1  models  (Figure  8a, b)  disappears  in 
Figure  8e.  which  is  the  also  the  case  of  the  D  =  r  models 
(Figure  8c, d) .  The  terminal  nitrogen  distribution  of  all 
the  atoms  of  the  D  —  r,  rt  =  °°  model  is  closest  to  that 
of  the  hydrated  model,  only  slightly  broadened. 

From  Table  2  and  Figures  6—8,  the  best  condition  to 
mimic  simulations  with  water  within  the  limitation  of 
our  hydration  model  is  D  —  r,  rt  =  °°  at  neutral  and  low 
pH  and  D-  1,  rc  =  llAat  high  pH.  Still,  the  D=  r,  rc 
=  °°  model  at  high  pH  is  reasonable  if  we  keep  in  mind 
that  the  overall  structure  is  more  condensed  (a  narrower 
and  left-shifted  profile  in  the  atom  density  distribution) 


than  the  ICN1  model.  The  atom  density  profile  of  the 
D=  r,  rc  =  °°  model  is  narrower,  but  the  density  profile 
of  terminal  nitrogen  is  very  close  to  corresponding 
results  from  simulation  with  water.  In  particular,  the 
terminal  groups,  which  would  be  essential  for  chemical 
and  biological  functionality,  show  similar  behavior 
which  validates  the  use  of  waterless  simulation  for  the 
current  purpose.  Therefore,  we  will  use  PAMAM  den- 
drimer  simulations  including  all  long-range  interactions 
with  distance-dependent  dielectric  constant  as  a  sub¬ 
stitute  for  the  explicit  water  simulation  for  the  remain¬ 
der  of  this  paper. 

We  studied  the  PAMAM  dendrimer  configuration  as 
a  function  of  pH  and  dendrimer  generation.  The  final 
configurations  of  PAMAM  dendrimer  generation  6,  after 
100  ps  of  MD  simulation,  at  different  pH  conditions,  are 
shown  in  Figure  9.  The  final  configuration  of  dendrimer 
at  high  pH  is  globular  and  loosely  compact  (Figure  9a). 
From  the  previous  simulation  results  and  discussion, 
we  know  the  overall  shape  is  probably  more  compressed 
than  simulations  with  water  molecules.  With  the  pro¬ 
tonation  of  the  primary  amines  at  neutral  pH,  the 
dendrimer  configuration  is  more  bulky,  bearing  large 
internal  cargo  spaces  due  to  electrostatic  repulsion 
(Figure  9b).  The  final  configuration  at  low  pH  presents 
a  higtily  ordered  structure  with  extensively  stretched 
dendrimer  branches,  generally  producing  an  open  struc¬ 
ture  (Figure  9c).  These  conformational  changes  by  pH 
can  be  applied  to  the  drug-loading  scheme  in  the 
dendrimer.  Analytic  comparison  is  done  by  the  atom 
density  distribution  and  the  density  of  terminal  nitrogen 
atoms  as  shown  in  parts  a  and  b  of  Figure  10,  respec¬ 
tively.  As  pH  changed  from  high  to  low  pH,  the  interior 
atom  density  is  reduced,  yielding  additional  interior 
space.  There  is  no  significant  initial  configuration 
dependency  of  the  ICN1  (blue)  and  ICN2  (red)  model 
(Figure  10). 

Dendrimer  conformation  change  was  investigated  as 
a  function  of  generation  number.  Atom  density  and  the 
density  of  terminal  nitrogen  atoms  of  the  PAMAM 
dendrimer  as  a  function  of  dendrimer  generation  at 
high,  neutral,  and  low  pH  are  presented  in  Figure  11. 
The  atom  density  is  normalized  by  total  atom  number 
of  G6  dendrimer,  and  tire  density  of  terminal  nitrogen 
atoms  is  normalized  by  the  total  number  of  terminal 
nitrogen  atoms  of  G6  dendrimer  at  each  pH  condition. 
As  the  generation  increases,  the  density  profile  moves 
to  larger  radial  distance  from  the  center  of  mass.  Table 
3  presents  maximum  radius  rmax  and  maximum  density 
radius  rd  of  all  these  graphs,  in  addition  to  Rc  values 
for  all  models. 

Figure  1  la  shows  atom  density  profile  over  radial 
distance  from  center  of  mass  of  genera tion-n  dendrimer 
at  high  pH.  The  density  graph  of  generation  n  is 
stacking  up  on  that  of  generation  n  —  1 .  Generation  n 
—  1  parts  are  conserved,  adding  the  additional  density 
in  the  outer  radial  region.  This  means  the  atom  distri¬ 
bution  at  high  pH  is  stable,  showing  that  the  dendrimer 
at  high  pH  preserves  compact  structure  as  far  as  branch 
stiffness  permits.  In  the  region  of  r  >  rd,  \dp/'dr\  (asymp¬ 
totic  value)  sharply  increases  from  generation  4.  Since 
the  density  profile  is  radial,  only  the  profile  at  r  >  rd 
can  be  compared  as  polymer  brush28  analogy.  Our 
graphs  present  similar  dropping  shape  with  the  results 
from  coarse-grain  MD  simulations  of  grafted  polymer 
brush:29  (1)  generations  2  and  3:  low  surface  density: 
(2)  generations  4,  5,  and  6:  high  surface  density. 
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Figure  9.  Final  configuration  of  G6  PAMAM  dendrimer  after  100  ps  MD  simulations  (a)  at  high  pH  (no  amine  is  protonated), 
(b)  at  neutral  pH  (primary  amines  are  protonated),  and  (c)  at  low  pH  (all  amines  are  protonated). 
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Figure  10.  (a)  Total  atom  density  and  (b)  the  density  of 
terminal  nitrogen  atoms  of  G6  PAMAM  dendrimer  at  different 
pH  conditions  from  ICN1  (blue  line)  and  ICN2  (red  line).  As 
pH  decreases,  the  density  profile  moves  to  a  larger  radial 
distance  from  the  center  of  mass. 

However,  |9p/i)ri  at  r  >  rd  (except  for  the  tail  region)  of 
G4  dendrimer  in  Figure  lib  (the  density  profile  of 
terminal  nitrogen  atoms)  is  between  those  of  the  G2  and 
G3  dendrimers  and  those  of  the  G5  and  G6  dendrimers. 
Therefore,  dendrimers  higher  than  generation  5  seem 
to  demonstrate  a  dense-packing  structure  relative  to  the 
intermediate  G4  dendrimer.  Interestingly,  the  terminal 
nitrogen  density  graph  of  generation  n  also  presents 
stacking  up  feature  on  that  of  generation  n  —  1.  This 
means  a  longer  branched  chain  end  has  the  same 
probability'  of  being  in  the  inner  region  of  the  dendrimer. 
The  branched  dendrimer  chains  should  be  flexible 
enough  for  this  back-folding  to  occur.  The  mean  shift  of 
rd  from  atom  density  profile  to  terminal  nitrogen  density 


one  for  all  generations  is  about  123%  at  high  pH.  These 
shifts  do  not  show  generation  dependency. 

The  density'  profile  at  neutral  pH  shows  a  different 
generation  dependency  as  shown  in  Figure  11c, d.  The 
atom  density  of  generation-/?  dendrimer  is  higher  than 
that  of  generation  n  —  1,  in  the  region  around  rd  of 
generation  n—  1.  Higher  generation  dendrimers  provide 
a  denser  intermediate  region  with  a  sparse  interior.  If 
a  dendrimer  is  used  for  drug  loading,  it  may  be  difficult 
to  put  the  drug  cargo  into  the  interior  of  higher 
generation  dendrimer  at  neutral  pH.  but  it  could  be  a 
good  carrier  once  drugs  are  loaded.  The  density  of 
terminal  nitrogen  atoms  appears  at  broad  radial  range 
for  each  generation  at  neutral  pH.  This  implies  that 
significant  branch  back-folding  occurred  at  neutral  pH, 
in  addition  to  the  major  peripheral  distribution  of  the 
terminal  groups.  The  density  profiles  of  terminal  nitro¬ 
gen  atoms  show  multiple  local  peaks  for  all  generations. 
There  seems  to  be  preferential  radii  for  branched 
terminal  atoms  to  locate.  |  3p/9/j  at  r  >  rd  increases 
gradually  as  generation  grows  in  both  total  atom  and 
terminal  nitrogen  atom  density  profiles.  The  mean  shift 
of  /'a  from  atom  density  profile  to  terminal  nitrogen 
density  one  is  about  1 1 4%  at  neutral  pH. 

Contrary  to  the  neutral  pH  case,  the  atom  density 
profile  of  generation-// dendrimer  at  low  pH  (Figure  1  ie) 
is  lower  than  that  of  generation  n  —  1,  in  the  whole 
region  of  generation  n  -  1 ,  As  generation  increases,  the 
dendrimer  molecule  becomes  less  dense  in  the  interior 
and  extends  out  from  the  center  of  mass.  Therefore, 
higher  generation  dendrimer  at  low  pH  may  be  the  best 
to  diffuse  in  and  out  small  noncharged  molecules  such 
as  drugs.  |9 p/3/j  at  r  >  rd  increases  evenly  as  generation 
grows  in  atom  density  profiles.  The  density  profiles  of 
terminal  nitrogen  atoms  are  localized  on  the  periphery 
of  dendrimer  molecules  as  shown  in  Figure  1  If.  At  low 
pH,  protonated  inner  amines  prevent  the  protonated 
terminal  amines  from  folding  back  into  the  inside  of  the 
dendrimer,  leaving  terminal  amines  on  the  periphery 
of  the  dendrimer.  The  minimum  radii  rmj„  for  nonzero 
terminal  density  are  16,  19.9,  25.9,  17,  and  25.4  A  for 
generations  2  to  6,  respectively.  The  ratios  of  (rd  -  rmin 
)/(.rmax  -  fij)  at  generations  2-6  are  1.2,  1.8,  1.6,  4.3, 
and  4.2.  These  kinds  of  G5  and  G6  dendrimer  graphs 
are  reminiscent  of  the  free  end  probability  graphs 
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Figure  11.  Density  profile  of  the  PAMAM  dendrimer  at  high  pH  (a.  b).  neutral  pH  (c.  d).  and  low  pH  (e,  f)  as  a  function  of 
dendrimer  generation:  generation  2  (black).  3  (pink),  4  (green),  5  (red),  and  6  (blue).  The  atom  densities  (left  panel:  a,  c,  e)  are 
normalized  by  total  atom  number  of  G6  dendrimer  and  the  densities  of  terminal  nitrogen  atoms  (right  panel:  b,  d.  f)  are  normalized 
by  the  number  of  toted  terminal  nitrogen  atoms  of  G6  dendrimer  at  each  pH.  As  the  generation  increases,  the  density  profile 
moves  to  a  larger  radial  distance  from  the  center  of  mass. 


Table  3.  Maximum  Radius  rlnax  and  Maximum  Density 
Radius  r,i  of  Figure  11,  in  Addition  to  /?c  Values  Tor  All 
Models 


Al  (A) 


generation 

pH 

fuiax 

(A) 

atom 

density 

density  of 
terminal  N 

Rc  (A) 

2 

low 

27.7 

19.2 

22.4 

16.6  ±0.1 

neutral 

25.3 

17.0 

19.6 

14.5  ±0.3 

high 

19 

9.3 

10.5 

8.4  ±0.1 

3 

low 

35.8 

26.3 

30.1 

22.8  ±  0.2 

neutral 

34.6 

22.7 

27.8 

19.7  ±0.1 

high 

21.6 

11.2 

14.3 

11.6  ±0.2 

4 

low 

44.7 

33.9 

37.6 

29.9  ±0.1 

neutral 

42.7 

31.0 

36.0 

26.7  ±0.2 

high 

27.4 

15.3 

19.2 

14.8  ±0.1 

5 

low 

54.1 

42.3 

47.1 

38.0  ±0.1 

neutral 

51.7 

39.4 

42.5 

32.8  ±  0.7 

high 

32.8 

19.5 

24.0 

18.3  ±0.3 

6 

low 

64.2 

52.3 

56.8 

46.8  ±  1.0 

neutral 

61.8 

47.3 

51.8 

41.3  ±0.1 

high 

44.1 

24.8 

30.7 

24.2  ±  0.4 

obtained  from  coarse-grain  MD  simulation  of  polymer 
brush  at  higher  surface  densities.29  |9p/9/1  at  r  >  of 
terminal  nitrogen  density  graphs  increases  gradually 
as  generation  grows  except  generations  2  and  3,  and 
the  mean  shift  of  rd  from  atom  density  profile  to 
terminal  nitrogen  density  is  about  1 12%. 


A  simple  model  for  polymer  brushes  gives  polymer 
height  h  from  free  energy  minimization  such  as28 

h  ~  N[wual)Vi  (3) 

where  a  is  the  Kuhn  length,  Nis  length  measured  by 
a,  w  is  the  excluded-volume  parameter,  and  o  is  the 
surface  grafting  density.  Our  model  of  PAMAM  den¬ 
drimer  is  spherical  with  branch  functionality  being 
two.14  It  is  difficult  to  properly  correlate  a  property  with 
o  for  comparison  purposes.  However,  Afcan  be  relatively 
correlated  to  the  value  of  generation  —  1.  As  shown  in 
Figure  12a.  EDA  core  part  has  roughly  half  of  the 
branch  length.  If  we  ignore  the  EDA  core,  the  radial 
difference  between  n  and  n  —  1  branch  points  can 
analogue  to  a  as  shown  in  Figure  12b;  then  the  radial 
distance  L  measured  by  radial  a  equals  generation  — 
1 .  From  Table  3,  we  could  obtain  the  rinax  and  L  relation 
for  each  pH  condition  such  as 

Aft  (4) 

The  values  of  exponent  a  are  0.97,  1.04,  and  0.99  at 
high,  neutral,  and  low  pH,  respectively.  Even  though 
the  dense  packing  seems  to  start  from  generation  5  at 
high  and  low  pH  on  the  basis  of  the  shape  of  terminal 


Macromolecules,  Vol.  35,  No.  11,  2002 


Figure  12.  (a)  Chemical  structure  of  GO  PAMAM  dendrimer 
with  EDA  core,  (b)  Schematic  drawing  of  G2  PAMAM  den¬ 
drimer  where  tiie  circles  represent  central  structure  (a). 
Denoted  numbers  are  the  radial  distance  L  (=  generation  - 

D- 


Figure  13.  Volume  change  between  neutral  and  low  pH  as  a 
function  of  generation.  Volumes  w'ere  calculated  as  Rq  assum¬ 
ing  a  spherical  dendrimer  model.  Lines  are  added  to  show 
linear  regions. 

nitrogen  distribution  and  literature,14  rmax  grows  quickly 
as  L  increases.  On  the  other  hand,  total  number  of 
branch  arms  in  a  dendrimer  can  be  defined  as  N  (28, 
60,  124,  252,  and  508  for  G2  to  G6  dendrimer,  respec¬ 
tively).  In  this  case,  the  exponents  are  0.29.  0.30,  and 
0.29  at  high,  neutral,  and  low  pH,  respectively.  Obvi¬ 
ously,  rmax  grows  very  slowly  as  A/increases.  which  gives 
the  small,  defined  size  of  dendrimers.  while  the  chain 
dimension  Rc  of  random-coiled  polymer  grows  as  A/1''2. 

Previous  gene  delivery  experiments  found  a  good 
amount  of  gene  expression  from  dendrimer  generations 
greater  than  5. 3  One  characteristic  of  a  dendrimer  as  a 
gene  delivery  vector  is  thought  to  be  the  ability  to  escape 
the  endosomal  compartment.  This  may  result  from 
further  protonation  of  the  dendrimer  at  endosomal  pH 
5,  as  compared  to  pH  7.4  outside  the  cell  and  in  the 
cytosol,  causing  endosomal  disruption  due  to  volume 
change  driven  by  osmotic  forces.30  The  volume  change 
between  neutral  and  low'  pH  is,  therefore,  of  interest. 
Figure  13  represents  the  volume  change  between  neu¬ 
tral  and  low  pH  at  each  generation.  Volumes  were 
calculated  as  Rc  assuming  a  spherical  dendrimer  model. 
A  transition  of  volume  change  happens  at  generation 
4.  This  feature  was  also  observed  from  Figure  lib. 
There  should  be  a  distinct  improvement  in  cytosolic 
delivery  from  generation  4  to  5.  Therefore,  we  can 
account  for  the  generation  dependency  of  gene  expres¬ 
sion  based  on  this  endosomal  disruption  hypothesis. 

Conclusion 

Extensive  MD  simulations  of  G2  PAMAM  dendrimer 
with  explicit  water  molecules  were  performed  from  two 
different  initial  configurations  over  a  range  of  pH 
conditions.  The  effects  of  long-range  interaction  param¬ 
eters  such  as  interaction  distance  and  dielectric  constant 
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were  investigated  in  simulations  of  the  same  dendrimer 
models  without  water  in  order  to  find  an  optimal  set  of 
conditions  for  waterless  simulation.  After  comparison 
of  results  for  hydrated  and  nonhydrated  G2  PAMAM 
dendrimers  over  a  range  of  simulation  parameter  set¬ 
tings,  we  found  the  optimal  condition  was  a  distance- 
dependent  dielectric  constant  without  a  nonbonded 
energy  cutoff.  The  general  simulation  results  of  the  G2 
PAMAM  dendrimers  with  these  parameters  corroborate 
the  simulations  done  with  10  A  explicit  water  layers, 
especially  in  the  charged  dendrimer  models  at  neutral 
and  low  pH.  At  high  pH,  the  overall  conformation  is  a 
slightly  compressed  under  these  long-range  interaction 
conditions  compared  to  the  hydrated  model. 

PAMAM  dendrimers  of  generations  3—6  were  con¬ 
tinuously  simulated  using  optimal  parameter  settings 
for  waterless  simulation  over  a  range  of  pH  conditions. 
The  radial  atom  density  of  each  dendrimer  was  deter¬ 
mined  for  each  pH  value.  Globular  and  loosely  compact 
structures  for  generation-/?  dendrimer  at  high  pH  (>  10) 
showed  a  radial  atom  density  profile  which  conserves 
the  generation  n  —  1  interior  while  adding  the  ad¬ 
ditional  outer  layer.  Highly  ordered  extended  genera¬ 
tion-;?  dendrimer  structures  at  low  pH  (<4)  produce  an 
atom  density  profile  that  is  reduced  in  the  region 
corresponding  to  generation  n  —  1.  This  hollowing  of 
the  interior  as  generation  increases,  within  our  simula¬ 
tion  range,  can  provide  increased  access  to  the  interior 
by  other  chemical  agents.  By  contrast,  significant 
branch  back-folding  occurring  at  neutral  pH,  in  addition 
to  the  major  peripheral  distribution  of  the  terminal 
groups,  yields  an  atom  density  profile  of  generation-/? 
dendrimer  that  is  higher  in  the  region  around  the 
maximum  density  of  the  generation  /?  -  1 .  This  suggests 
that  the  higher  generation  dendrimer  provide  a  cavity 
surrounded  by  dense  atom  populations,  which  should 
be  a  more  stable  agent  carrier.  The  volume  differences 
between  neutral  and  low  pH  calculated  from  Rc  showed 
a  significant  increase  beginning  with  generation  5. 
These  structural  differences  are  possibly  responsible  for 
the  improvement  in  gene  expression  at  G5  dendrimer 
as  a  gene  delivery  system. 
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